Bootstrap tomography of high-precision pulses for quantum control 
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Long-time dynamical decoupling and quantum control of qubits require high-precision control 
pulses. Full characterization (quantum tomography) of imperfect pulses presents a bootstrap prob- 
lem: tomography requires initial states of a qubit which can not be prepared without imperfect 
pulses. We present a protocol for pulse error analysis, specifically tailored for a wide range of the 
single solid-state electron spins. Using a single electron spin of a nitrogen-vacancy (NV) center in 
diamond, we experimentally verify the correctness of the protocol, and demonstrate its usefulness 
for quantum control tasks. 
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Coherent manipulation of single and few electron spins 
has recently been achieved in several solid-state sys- 
tems such as quantum dots and diamond defect cen- 
ters. Such systems are promising candidates for quan- 
tum information processing P, 0] , precise metrology [3] 
and ultra-sensitive magnetometry They also present 
an excellent testbed for studying the fundamental prob- 
lems of quantum dynamics of open systems [HI-ItII. High- 
speed manipulation of the system's quantum state can 
be achieved by using microwave or optical pulses 
which must be fine-tuned to provide a high degree of fi- 
delity. For example, sequences of quantum control pulses 
can be used to exte nd the coherence time via dynamical 
decoupling 



11H14, 



16[. For long sequences, even small 
errors in the pulses will destroy the coherence that one 
attempts to preserve IJ, |15| and may even lead to arti- 



ficial saturation 17, 18|. Therefore, precise characteriza- 
tion of errors is essential for successful implementation of 
complex quantum control protocols. With known errors, 
composite pulses and/or special pulse sequences can be 
chosen to mitigate the problem. 

Complete information on the action of a pulse can in 
principle be gained with quantum process tomography 
(QPT) . However, QPT of an imperfect pulse re- 
quires preparation and measurement of a complete set of 
reference states, whereas in many solid-state qubit sys- 
tems (e.g. quantum dots, diamond defect centers, super- 
conducting circuits) only one state can be prepared re- 
liably (without the imperfect pulses), and only one ob- 
servable can be directly measured. All other states can 
be prepared only with the imperfect pulses themselves, 
and therefore have errors [2^ . This presents a bootstrap 
problem: the reference states contain the very same er- 
rors that we want to determine. 

The problem of pulse error analysis has been studied 
extensively in the areas of NMR and ESR |2l| - |24| . How- 
ever, single electron spins in solid-state settings present 
new opportunities and challenges, and call for new ap- 
proaches tailored at the specific demands of these sys- 



tems. The driving pulse field can be tightly confined 
in the vicinity of the target spin. The resulting strong, 
nanosccond-timescale pulses enable fast spin manipula- 
tion, but the standard pulse error analysis 2ll424l| used 
in NMR becomes inapplicable. At strong driving, the 
spin dynamics changes noticeably [lo| . The non-secular 
terms in the rotating frame can become important. The 
ac-Stark and Bloch-Siegert shifts can significantly detune 
the pulse frequency from resonance 10] and tilt the ro- 
tation axis towards the z-axis. Also, the pulse edges con- 
stitute a much larger fraction of the short pulse, and the 
driving field at the edges varies much faster and stronger 
than in typical NMR pulses. The resulting errors [10| 
(e.g. tilting of the rotation axis) can go beyond the stan- 
dard treatment [25], and can not always be removed by 
symmetrizing the pulse shape. 

Also, typical NMR systems have long coherence times 
that exceed the pulse width by orders of magnitude. The 
standard tune-up protocols 2l|-|2j] exploit this advan- 
tage, and use sequences with tens or hundreds of pulses to 
achieve outstanding precision in pulse parameters. But 
single solid-state electron spins are dephased faster, on 
a timescale of microseconds down to tens of nanosec- 
onds Q . After only tens of pulses the signal becomes a 
complex mixture of pulse errors and decoherence 17. 18|. 
To ensure a reliable measurement of the errors, the pro- 
tocol for single electron spins must be short, a few pulses 
at most. 

Here, we present a systematic approach to pulse char- 
acterization for single solid-state electron spins, which 
is usable at shorter coherence times and much stronger 
driving power compared to traditional NMR systems. 
The proposed protocol contains four series of measure- 
ments, each having only 1-3 pulses, thus minimizing the 
effect of decoherence. The measured signal quantifying 
the pulse errors grows linearly with the errors to ensure 
a good accuracy for small errors. Also, the signal is zero 
for zero errors for good relative accuracy. The protocol 
determines all pulse errors: the rotation angle and all 
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three components of the rotation axis [25|- We experi- 
mentally demonstrate the protocol on a single spin of a 
nitrogen- vacancy (NV) defect center in diamond. By de- 
liberately introducing known pulse errors, we verify the 
accuracy and self-consistency of the protocol, and use it 
to significantly increase the fidelity of QPT. 

Our goal is to determine the parameters of four pulses, 
TTx, TTy, 7r/2x, and 7r/2y applied to a two-level quantum 
system (tt^ denotes a rotation by an angle tt around the 
X-axis in the rotating frame; other notations are anal- 
ogous). These pulses allow implementation of univer- 
sal decoupling XY sequences ll|, [l^l , full tomography of 
the density matrix, and universal single-qubit gates [l9|. 
We assume that the pulse errors are reasonably small, 
and consider only the first-order terms in these quanti- 
ties (since we want the signal to grow proportionally to 
errors). We also assume that the pulse width tp is small 
in comparison with the dephasing time Tj*; in this case 
the impact of decoherence is of second order, (ip/Tj*)^, 
and is negligible for short sequences Under this as- 
sumption the evolution of a spin during the pulse can be 
described a unitary rotation. For example, for S = 1/2, 
the evolution (in the rotating frame) during an imperfect 
TTx pulse is given by 



Ux 



-4) ~ i{ax + eyCTy + ezO-z), (1) 



where (Jx,y,z are the Pauli matrices, the rotation angle 
error is 2(j) and the rotation axis n has small components 
Uy — Cy and Uz — tz- Similarly, a. 'n:/2x pulse U'^ has 
the rotation angle error 2(j)' , and the small rotation axis 
components e'y and e'^ along y and z, respectively. Note 
that in general two 7r/2 pulses do not yield the same 
evolution as one tt pulse due to errors introduced by the 
pulse edges. Analogous parameters for y-pulses will be 
denoted as 2%, v^, and Vz (angle and axis errors for Try), 
and 2x', v'^, and v'^ (angle and axis errors for it/2y)- 

The bootstrap protocol shares ideas with the standard 
QPT, and with the NMR tune-up sequences. Before each 
measurement, the spin is in the state | t)i and the mea- 
sured signal is z\tp) , where 1-0) is the wavefunction 
after the pulse. An imperfect pulse Uj can be repre- 
sented as a product Uj = U^^^Vj ~ uf\l — iKj), where 
AO) 



7fV,«[/f (1- 

Uj"^' is a corresponding ideal rotation and the Hermitian 
operator Kj is proportional to small pulse errors. Apply- 
ing two pulses Ui and U2 in succession, we obtain up to 
linear order in 



U21 = U2U1 « C/i°^C/|"^ - lU^Ki - iK2Ul 



(0), 



r(0) 



(2) 



and the terms ui^^ Ki and K2Uf''' contain different ma- 
trix elements of the operators Ki and K2. E.g., if Ui 
and U2 are the (imperfect) 'k/2y and 'n/2x rotations, the 
signal detected after this sequence, S21 = Tr((T2[/2i| t)(t 
|?7|j^), contains a linear combination of the matrix ele- 
ments (t \Ki\Y) and (t \K2\X) (where \Y) = \^)+i\\.) 



TABLE I: Summary of the bootstrap protocol: pulse se- 
quences (read from right to left) and the resulting signals ex- 
pressed in terms of the error parameters. Blocks of sequences 
are separated by horizontal lines. 



Sequence 


Signal 


7r/2x 


-20' 


7r/2y 


-2x' 


■k/2x-i^x 


2(0 + 0') 




2(X + X') 


■RY--K /2x 


-2v^ + 20' 


7rx-vr/2y 


+ 2x' 


■k/2y-it/2x 


-4 - 4 - 4 - 


7r/2x-7r/2y 




7r/2x-7rx-7r/2y 


-Ej/ + £z + — + 2<iy 


7r/2y-7rx-7r/2x 


-ey ~ ez+v^+Vz+ 2ey 


TV I2x-'^Y-'^ I2y 


e'y ~ e'z - v'^ + v'z + 2v^ 


7r/2y-7ry-7r/2x 


e'y + - v'^ - + 2Vj: 



and \X) = I '1^) -I- I D). Combining different pulses, we 
obtain a sufficient number of such linear combinations of 
various matrix elements of Kj to uniquely determine all 
of them. A general approach to bootstrap tomography 
can be formulated in the language of QPT, by expanding 
the operation element operators [l^ in terms of small er- 
rors, and various bootstrap protocols applicable to more 
complex systems can be designed in a similar manner. 
Here, we focus on a protocol for a single two-level sys- 
tem. 

The proposed protocol is summarized in Table HI It 
consists of three blocks of measurement sequences. For 
each sequence the measured signal is given in terms of 
the error parameters. The first block, with two single- 
pulse sequences, yields the rotation angle errors for the 
7r/2 pulses. This information is then used in the second 
block, consisting of four two-pulse sequences, to find the 
rotation angle errors and the components of the rotation 
axis along z for the tt pulses. The third block has six 
multi-pulse sequences, yielding six signals that are lin- 
early related to the remaining six pulse error parameters. 
This linear system is underdetermined, since the whole 
system of pulses is invariant under rotations around the 
z-axis. We may put e'y = 0, taking the phase of the 7r/2x 
pulse as the x direction in the rotating frame. This fixes 
all other directions, and all errors are uniquely deter- 
mined. No unphysical results, typical for experimental 
implementations of standard QPT [2^ (2^ , appear in this 
bootstrap protocol. 

We now demonstrate and verify the protocol experi- 
mentally by applying it to a single solid-state spin sys- 
tem. We use the spin of a single Nitrogen- Vacancy (NV) 
center, which is a defect in diamond composed of a sub- 
stitutional nitrogen atom with an adjacent vacancy {26j. 
The NV center's spin can be optically initialized and read 
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FIG. 1: (Color online) Experimental verification of the boot- 
strap protocol by introducing varying pulse errors. Duration 
of the 7r/2-pulses (7r-pulses) is 5 ns (9 ns). (a) Measured er- 
ror parameters for diflterent phases $ of the 7r/2y-pulse. The 
frequency of the driving field is set at 2.4605 GHz. (b) Mea- 
sured error parameters for various frequencies of the driving 
field. Error bars everywhere are smaller than the symbol size. 



out [26j . The experiments are performed in a home- 
built confocal microscope at room temperature. NV 
centers in nanocrystals are prepared on a chip with a 
UthographicaUy-defined waveguide ahowin g fa st and pre- 
cise spin rotations by magnetic resonance [25| . 

We controUably introduce two types of pulse errors, 
and use the bootstrap protocol to extract their values. 
In the first experiment, we vary the phase $ of the nom- 
inal 7r/2y-pulse between —30° and 30° from its nominal 
value. This way, we are changing the error parameter 
v'^ = sin<i> w $(rad) while leaving all other error param- 
eters constant. Figure [ija) shows the results from the 
bootstrap protocol that clearly support this expectation. 

In the second experiment we detune the microwave 
excitation away from the qubit transition frequency, 
thereby varying the z-components of the rotation axis 
for all pulses. As shown in Fig. [Ijb), the extracted error 
parameters Vz, v'^, Cz, and strongly change (roughly 
linearly) with the detuning as expected, while the other 
error parameters stay virtually constant. If these rota- 
tion axis errors were arising only from the bulk of the 
pulse, they would all show the same dependence on the 
detuning, and all four curves would have the same slope. 
Instead we observe that the errors of the nominal 7r/2- 
pulses vary about twice as much as the errors of the nom- 
inal TT- pulses. This indicates that the errors originate 
largely from the pulse edges. Since the edges are the 
same for all pulses, they have a relatively larger impact 
on shorter pulses. The data in Fig.[lja)-(b) demonstrate 
that the bootstrap protocol is indeed an efficient and re- 
liable tool for extracting pulse errors. 

Due to experimental limitations it may be impossible 
to cancel all errors at once. In that case, the choice of 
the optimal working point involves a trade-off, and pre- 
cise knowledge of the pulse errors becomes particularly 
important. For example, when performing QPT, a set 
of the reference states is prepared using the pulses ttx, 



7r/2x, and 7r/2y. These states are acted upon by the 
process, and rotated to the readout basis before mea- 
surement [1^. The operation elements of the quantumn 
process are expanded in the basis Eq = I, Ei = ax, 
E2 — cTy, and i?3 — Uz, and the process is completely 
characterized by the 4x4 expansion matrix x • When 
systematic pulse errors are present, the prepared initial 
states differ from the reference states, and the read-out 
is also performed in the incorrect basis, yielding an in- 
correct matrix x- But with pulse errors known, the raw 
measured data can be transformed into the correct basis 
prior to the standard QPT data processing [3 US [23. 

As a demonstration, and as a check of self-consistency 
of the bootstrap protocol, we perform QPT while intro- 
ducing the same pulse errors as in Fig. 1. We show that 
with the pulse errors deduced with the protocol, the QPT 
results can be corrected. The comparison between raw 
and corrected data below is designed to use no a priori 
assumptions about correctness of the bootstrap protocol. 

First, we take the (imperfect) tty pulse as an example 
of a quantum process. We introduce errors in the QPT 
procedure by changing the phase $ of the nominal 7r/2y- 
pulse from —30° to 30°. We first determine the reference 
process which corresponds to $ = 0; the corresponding 
experimental settings and the pulse error parameters are 
marked in Figs. 1 and 2 by dotted lines. We perform 
QPT on this reference process, and the resulting refer- 
ence matrix xo is calculated in two ways: (i) using the 
raw uncorrected data, i.e. assuming that the pulses used 
for QPT are ideal (we denote this matrix as Xo)i and (ii) 
using the data corrected for the known pulse imperfec- 
tions (the resulting matrix is Xo)- Next, we vary $, and 
use artificially deteriorated 7r/2y pulses to determine the 
matrix x of the quantum process. This matrix is also 
determined in two way, by using raw experimental data 
(matrix x^), and by correcting the data for the known 
pulse errors (matrix x^)- For each value of we com- 
pare the raw-data matrices x'' and Xo on 0116 hand, and 
the corrected matrices x'^ and Xo on the other. 

The process we are studying does not depend on the 
phase of the nominal it /2y pulse. Thus, ideally, the ma- 
trices Xo and x should be the same. To quantify the dif- 
ference between xo and X) we use two distance measures. 
One is the process fidelity [3| F = Tr[xox]i which de- 
pends quadratically on the pulse errors. The other mea- 
sure is the Hilbert-Schmidt 2-norm pf II2 = ^/T^\mW] 
of the difference matrix M = x — Xo- This norm is linear 
in, and thus more sensitive to, the pulse errors [2^. 

In Figs. [2][a)-(b), orange squares show the values of 
F and ||M||2 for the corrected-data matrices Xo and x'^- 
The expectation that xo and x should coincide is con- 
firmed with excellent precision. Almost independently of 
<&, the fidelity remains above 99%, and | \M\ [2 stays small. 
This is not so for the raw-data matrices Xo and x' (blue 
squares). The neglected phase error of the nominal 7r/2y 
pulse makes the matrix x"^ inaccurate, so F and ||-A/||2 
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FIG. 2: (Color online) Correction of pulse errors in Quantum 
Process Tomography using the bootstrap protocol, (a) Fi- 
delity (F) and (b) the 2- norm distance ||M||2 between the 
process measured at finite introduced 7r/2y phase error and 
the process matrix measured at zero introduced error. The 
process is a Try-pulse with zero introduced error, (c) F and 
(d) ||M||2 between the measured process and the actual pro- 
cess (identity). All measures are calculated both for the un- 
corrected and for the corrected data. 



depend on $, with fidelity dropping by 8% for $ = 30°. 

In a second experiment, we detune the microwave ex- 
citation frequency away from the qubit transition, intro- 
ducing the errors e^, e'^, Vz, and v'^ (see Fig. [T}d) into all 
pulses. The tomographed process is the identity process, 
which is independent of the pulse carrier frequency. We 
perform QPT on this process for the same range of detun- 
ings as in Fig.[IJb). As above, we determine the corrected 
and the uncorrected matrices x'', and x'^, while the ma- 
trix xoi describing the identity process, is independent of 
the pulse errors. The results are shown in Fig. [irc)-(d). 
Again, the fidelities are high for the corrected data for the 
full range of introduced errors, while for the uncorrected 
data the fidelity has dropped by as much as 10%. The 
same behavior is seen for ||Af||2. Thus, even the effects 
of complex pulse errors introduced by detuning the fre- 
quency can be efhciently corrected using the information 
from the bootstrap protocol. 

Summarizing, we have developed an effective pulse er- 
ror analysis protocol tailored to the specific requirements 
of single solid-state spins. We have experimentally im- 
plemented and verified the protocol on a single electron 
spin in diamond. We have shown that the distortions of 
the tomography results, arising due to the pulse errors, 
can be corrected with the knowledge obtained from the 
bootstrap protocol. The methods described in this paper 
are applicable to many systems. They may help in accu- 
rate determination of the properties of different quantum 
processes, a key feature for the fields of quantum infor- 
mation processing, quantum metrology and fundamental 



studies of quantum decoherence. 
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